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Abstract 

CD In this paper, we consider the motion of the hydrodynamic Chaplygin sleigh, a planar rigid body 

in a potential flow with circulation around the body, subject to a nonholonomic constraint modeling 
a fin or keel attached to the body. We show that the motion of this system can be described by 
Euler-Poincare-Suslov equations on a central extension of the special Euclidian group SE(2), where 
the cocycle used to construct the extension encodes the effects of circulation upon the body. In the 
second part of the paper, we then discuss nonholonomic systems on central extensions of Lie groups. 
Oh where both the Lagrangian and the nonholonomic constraints are left invariant. We show that 

there is a one-to-one correspondence between invariant measures on the original group and on the 
extended group, and we use this result to characterize the existence of an invariant measure for the 
hydrodynamic Chaplygin sleigh. We finish with a qualitative discussion of the reduced dynamics. 
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1 Introduction and outline 

In this paper, we consider the motion of a planar rigid body surrounded by an irrotational perfect 
fluid with a given amount of circulation around the body. In addition, we assume that the body is 
equipped with a very effective keel or fin. To a good approximation, this fin imposes restrictions on the 
possible velocities of the body, and can therefore modeled by a nonholonomic constraint. In this way 
we obtain a low-dimensional model of an underwater vehicle with a rigid fin, which we have termed 
the hydrodynamic Chaplygin sleigh with circulation, after [lOJ where the corresponding model without 
circulation was described. 

History and overview of fluid-structure interactions. The motion of a rigid body in a potential 
fluid in the absence of external forces was first described by Kirchhoff |17] . His crucial observation 
was that the effect of the fluid on the body could be described entirely in terms of the added mass 
and added inertia terms, which depend on the geometry of the body only, and can be calculated 
analytically for a wide class of body shapes. Kirchhoff's solution was extended to the case of rigid 
bodies moving in potential flow with circulation by, among others, Chaplygin [6j and Lamb |19| . who 
derived the equations of motion for this system, provided an explicit integration in terms of elliptic 
functions, and described qualitative features of the dynamics. In recent years, these ideas have been 
extended to the case of rigid bodies interacting with point vortices [31^ [1] , vortex rings [32] and other 
vortical structures, and they have been used to describe underwater vehicles pO] and the motion of 
bio-organisms [HE]. A comprehensive overview of the history of these equations can be found in 
[2]. Fluid-structure models involving nonholonomic constraints have been considered in the aerospace 
engineering community [2H] , while robotic models for underwater vehicles taking into account the effects 
of circulation were considered in [15] . 

The hydrodynamic Chaplygin sleigh. The model described in this paper incorporates nonholo- 
nomic constraints into the Chaplygin-Lamb system described above. We refer to this model as the 
hydrodynamic Chaplygin sleigh, since in the absence of the fluid the nonholonomic constraint models 
the effect of a sharp blade in the classical Chaplygin sleigh problem |7j which prevents the sleigh from 
moving in the lateral direction. It is an interesting historic coincidence that the name of Chaplygin is 
linked both to the development of the Chaplygin-Lamb equations [6j as well as to the nonholonomic 
Chaplygin sleigh [7]. 

Consider a body frame {Ei, E2} where Ei is aligned with the fin of the body (see Figure [l|. The 



2 



Center of mass 

/_ 















Blade -^^^^ 


\ B 





Contact point 



Figure 1: The hydrodynamic Chaplygin sleigh consists of a rigid body moving in a potential fluid 
equipped with a very effective fin or blade, which prevents motion in the direction lateral to the 
fin. The body frame {Ei,E2} is taken so that Ei is aligned with the fin. The fin does not have 
to be aligned with the axes of the body, and we make no assumptions about the relative location 
of the center of mass, the geometric center of the body, and the contact point of the fin. We have 
depicted an elliptical body shape, but our results are valid for arbitrary convex shapes. 



equations of motion for the hydrodynamic Chaplygin sleigh with circulation have the following form: 

k = V2P1 — V1P2 — p{avi + PV2), 
pi = ujp2 — Kpv2 + pauj, (1.1) 
P2 = -ojpi + Kpvi + p(3uj + A. 

Here, p is the density of the fluid, (w, ^1,^2) are the angular and linear velocities of the body, {k,pi,p2) 
represent the corresponding impulses defined via the added inertia tensor, k and a, /3 are constants 
related to circulation, and A is a Lagrange multiplier determined uniquely from the nonholonomic 
constraint 

V2 = 0, (1.2) 

which precludes motion in the direction perpendicular to the blade. We emphasize that these equations 
are valid for arbitrary body shapes; the geometry of the body is encoded in the added inertia matrix of 



the body (see (2.10) below). 

If the constraint is not enforced and its Lagrange multiplier is set to zero, the corresponding equa- 
tions become those of Chaplygin and Lamb |19] . If instead the circulation vanishes (k = 0, and as a 



result, a = /3 = 0), the resulting system is the hydrodynamic Chaplygin sleigh treated in [TOi 



Hamiltonian structure of the Chaplygin-Lamb equations. The equations (1.1) have several 
interesting features. In section [3] we begin by focusing on the Chaplygin-Lamb equations: in other 
words, we consider a rigid planar body of arbitrary shape immersed in a potential flow with circulation. 



but we do not yet take the nonholonomic constraint (1.2) into account. In section 3.2 we show that 
the Chaplygin-Lamb equations are Hamiltonian with respect to the following Poisson bracket, defined 
on functions F,K of the angular and linear momenta p = {k,pi,p2)- 





{F,K},,^,f,{p) = {V^FY I p2 (V^K) - -p/3 • ((V^F) x (V^K)) . (1.3) 
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The first term in this expression is the Lie-Poisson bracket on the dual space se(2)*, while the second 
term is non-canonical, and is due to the presence of circulation. 



In section 3.3, we make a digression to give a geometric interpretation of the constants a, /3, and 
in particular their dependence on the circulation k. While these constants were introduced already 
by Lamb [19], their role remained somewhat mysterious. We introduce the cyclic centroid of a rigid 



body with circulation, and we show that a and /3 can be made zero by attaching the body frame to 



the cyclic centroid, leading to a corresponding simplification in the equations (1.1). Geometrically, the 
cyclic centroid is the center of the circle which best approximates the circumference of the rigid body, 
in the sense of complex approximation theory [9j (see Figure [4]) . It can then be shown that the cyclic 
centroid is the root (in the complex plane) of the first Faber polynomial of the rigid body, while the 
radius of the best-fitting circle is given by the conformal capacity of the rigid body. 

Circulation in terms of central extensions. In section |4] we show that the noncanonical Poisson 



bracket (1.3), defined on the dual of the Lie algebra of the Lie group SE(2), can be made into a canonical 
Lie-Poisson bracket by considering a central extension of SE(2) by M'^. The M^-valued cocycle used to 
construct this central extension is shown to consist of two parts: the first is an R- valued cocycle which 
models the effect of non-zero circulation k. This cocycle was introduced in [35j and gives rise to the 
oscillator group, a non-trivial central extension of SE(2) by M. The second cocycle takes values in M? 
and models the effect of the coefficients q, /3 in the Chaplygin-Lamb equations. In contrast to the first 
cocycle, this cocycle is exact, and we show that it can be "gauged away" by attaching the body frame 
of reference to the cyclic centroid. When nonholonomic constraints are present, this simplification is 
often undesirable, as it would lead to more complicated expressions for the constraints, so that we have 
to include the effect of nonzero a, /3. 

Nonholonomic mechanics on central extensions. In section [5] we incorporate nonholonomic 
constraints into our model. We begin by describing the class of LL systems on central extensions, 
nonholonomic mechanical systems for which the configuration space is the central extension of a Lie 
group, and where both the Lagrangian and the constraints are invariant under the left action of the 



extended group on itself. In the case of the equations (1.1), the underlying Lie group is the central 
extension of SE(2) by described previously, and the LL invariance expresses the invariance of the 



equations (1.1) and the constraint (1.2) under translations and rotations of the rigid body and the 
surrounding fluid. 

Given a Lagrangian and a set of nonholonomic constraints on a central extension G of Lie group 
G, we derive the corresponding reduced equations of motion, referred to as the Euler-Poincare-Suslov 
equations. We then turn to the question of whether these equations have an invariant measure. This 
problem has been considered for general EPS equations on Lie algebras of compact groups by Kozlov 
|18j . by Jovanovic |13| in the non-compact case, and by Zenkov and Bloch [36] for systems with 
nontrivial shape space. Using the criterion of Jovanovic, we show that the EPS equations on a central 
extension G admit an invariant measure if and only if the corresponding system on the original group 
G has an invariant measure, and we use this result to characterize the existence of invariant measures 



for the equations ( 1.1 ) in terms of the added moments of inertia of the hydrodynamic Chaplygin sleigh. 
In particular, we obtain that the existence of an invariant measure is independent of the value of the 
circulation k and the constants a, /?. 
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Qualitative features of the dynamics. We finish the paper with an overview of some features of 
the reduced dynamics. In particular, we show that there exists a one-parameter family of equilibria in 
the reduced phase space, and that the other solutions are asymptotic motions between the equilibria. 
A more detailed description of the reduced dynamics, including analytic expressions for the solutions, 
will be given in 



2 Kinematics of rigid bodies and ideal fluids 

We adopt Euler's approach to the study of the rigid body dynamics and consider an orthonormal body 
frame {Ei,E2} that is attached to the body. This frame is related to a fixed space frame {ei,e2} 
by a rotation by an angle 9 that specifies the orientation of the two dimensional body at each time. 
We will denote by x = (x, y) G M? the spatial coordinates of the origin of the body frame (see Figure 
[2]). The configuration of the body at any time is completely determined by the element g of the two 
dimensional Euclidean group SE(2) given by 

— sin X \ 

G SE(2). 




Center of mass 



Center of mass 




Origin of body frame 




Origin of body frame 



(a) Body frame is aligned with axes 
of symmetry of the body. 



(b) Arbitrary position and orienta- 
tion of the body frame. 



Figure 2: Two different choices of the body frame for an elliptical two-dimensional rigid body. In 
both cases the origin of the body frame does not coincide with the center of mass. 

We will often denote the above element in g £ SE(2) hy g = {Rg,x), where Re £ S0(2) is the 
rotation matrix determined by the angle 6. Let {vi, V2) € be the linear velocity of the origin of the 
body frame written in the body coordinates, and denote hy uj = 9 the body's angular velocity. They 
define the element in the Lie algebra se(2) given by 






—OJ 


Vl 


UJ 

















G sc(2) 



(2.4) 
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Explicitly we have 

6 = uj, vi = xcosO + ys'mO, V2 = —xsinO + y cosO. (2-5) 

For convenience, we will sometimes identify se(2) with as vector spaces, and denote ^ G se(2) as the 
column vector {uj, vi,V2)'^ G M^. The Lie algebra commutator takes the form 

[{UJ,VI,V2) , {UJ',V[,V2)]^^^^^ = (0, V2UJ' - U^v'^ , UJv\ - V^J ) ■ 

The kinetic energy of the body is given by 

= - ((Xg + mc? + mb'^)uj'^ + m{vf + w|) — mbujvi + maujV2) , 

where m is the mass of the body, (a, h) are body coordinates of the center of mass (see Figure [2]), and 
Xg is the moment of inertia of the body about the center of mass. It is a positive definite quadratic 
form on se(2) whose matrix is the body inertia tensor 

Xg + m[a? + 6^) —mb ma 
Ib = I -mb m I . (2.6) 

ma m 

For future reference, we give an explicit description of the dual space se(2)*. Since 5e(2) is isomorphic 
to M'^ and using the Euclidian inner product, we have that se(2)* = M^. A typical element /i is 
represented as a row vector /i = {k,pi,p2)- The duality pairing between fi and an element ^ = 
{u,vi,V2)'^ of se(2) is given by 

{fi, S,) = IJ.(, = kuj + pivi + P2V2. 

Lastly, we also mention that SE(2) acts on se(2) and sc(2)* by the adjoint and coadjoint actions, 
respectively. For the present paper, we will often need the infinitesimal version of the coadjoint action. 
This action is denoted by ad| : se(2)* — )• se(2)* for all ^ € se(2) and is given explicitly by 

ad|(/i) = {piV2 - p2Vi,ujp2,-ujpi), (2.7) 

where ^ = (w, vi, ^2)'^ and fi = {k,pi,p2). 

The dual space se(2)* is equipped with the (minus) Lie-Poisson bracket, which is given by 

for all functions F, K on 5e(2)*. In coordinates, we have 

-P2 Pi ' 

{^'^};e{2).(/^) = (Vm^)' I P2 I {V,K), 

0, 

where V^F is the gradient of F with respect to the variables {k,pi,p2) = /i. 
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The fluid flow at a given instant. Consider now the motion of the fluid that surrounds the body. 
Suppose that at a given instant the body occupies a region i3 C M^. The flow is assumed to take place 
in the connected unbounded region U := \ that is not occupied by the body. We assume that the 
flow is potential so the Eulerian velocity of the fluid u can be written as u = V^" for a fluid potential 
<I> : — )• M. Incompressibility of the fluid implies that <I> is harmonic, 

V^^ = on U. 



The boundary conditions for <I> come from the following considerations. On the one hand it is 
assumed that, up to a purely circulatory flow around the body, the motion of the fluid is solely due to 
the motion of the body. This assumption requires the fluid velocity V$ to vanish at infinity. Secondly, 
to avoid cavitation or penetration of the fluid into the body, we require the normal component of the 
fluid velocity at a material point p on the boundary of B to agree with the normal component of the 
velocity of p. Suppose that the vector {X, Y) G M? gives body coordinates for p. The latter boundary 
condition is expressed as 

TT- = ivi- ujY)ni + {v2+ujX)n2-, 

where n = (ni,n2) is the outward unit normal vector to B at p written in body coordinates. These 
conditions determine the flow of the fluid up to a purely circulatory flow around the body that would 
persist if the body is brought to rest. The latter is specified by the value of the circulation k around 
the body as we now discuss. 

The potential $ that satisfies the above boundary value problem can be written in terms of the 
body's velocities vi,V2,u}, in Kirchhoff form: 



^ = f 101 + V2(l)2 +UJX + ' 



(2. 



where (pi, i = 0,1,2, and x harmonic functions on U whose gradients vanish at infinity and 
satisfy: 



dn 



1,2, 



dB 



dx 
dn 



Xn2 - Yni, 



dB 



dn 



0. 



dB 



The potential (pQ is multi-valued and defines the circulatory flow around the body. The circulation k 
of the fluid around the body satisfies 



K = <7) u dl= f V4>o ■ dl, (2.9) 
'dB JdB 

and remains constant during the motion. Figure [3] shows the streamline pattern of the flow determined 
by the motion of an elliptical body for different values of k. 

Disregarding the circulatory motion, the kinetic energy of the fluid is given by 



|V($-</>o)||'dA 



u 



where dA is the area element in M? and p is the (constant) fluid density. We have subtracted the 
circulatory part from the velocity potential, as it would give rise to an infinite contribution to the fiuid 
kinetic energy [30j . 
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(a) K = 



(b) K = 27r 



(c) K = 10 



Figure 3: Stream line pattern for an ellipse moving on the plane for different values of the 
circulation k. The major and minor semi-axes of the ellipse are A = 2, B = 1. The body frame is 
aligned with the principal axes of the ellipse and the velocity of the body satisfies vi — V2 = 1 and 
w = 5. 



By substituting (2.8) into the above, one can express Tj- as the quadratic form 



Tt=\\Y. -^^^^^^j- + 2 ^^^^^ + ^^^^ I ' (2.10) 
\i,j=i i=i J 

where A^^, /Cjp, i, j = 1,2, and Xjr are certain constants that only depend on the body shape. Explicitly 
one has (see [19] for details), 

M'Jr = -p[ 4>i^dl = -p[ (pj^dl, i,j = 1,2, Xjr = -p[ x^dl 
JdB (J'"' JdB O*?! Jqi3 on 



jCjr = -p[ (j)i^dl = -p[ x^dl, i = 1,2. 



These constants are referred to as added masses and are conveniently written in 3 x 3 matrix form to 
define the (symmetric) added inertia tensor: 



that defines Tjr ctS cl quadratic form on se(2). 



2> /Cjr 



Example. For an elliptic rigid body with semi- axes of length A > B > 0, the added masses and 
moments of inertia take on a particularly convenient form. The kinetic energy of the fluid is given by 
(see [H]) 

T^ = ^(B^vl + A^vj+^-^^^^.^ 



2 \ ' ^ 4 
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where we have ignored the circulatory motion around the body. The corresponding added inertia tensor 
is thus given by 



Ijr = piT 




(2.11) 



We emphasize that this particular form of the added inertia tensor was derived under the assumption 
that the body frame is aligned with the symmetry axes of the ellipse {u = 0, r = s = in Figure [T]). 
When this is not the case, the added mass tensor is more complicated, and in particular need not be 



diagonal, as is shown in (5.40). 



3 Hamiltonian dynamics of rigid bodies in potential flow 

In this section, we describe the Kirchhoff equations for the dynamics of a rigid body in potential flow, 
and we review the Chaplygin-Lamb equations dealing with rigid bodies in the presence of circulation. 
We show how, in the absence of circulation, the equations of motion can be viewed as Lie-Poisson 
equations on the space se(2)*, the dual of the Lie algebra of the Lie group SE{2) of translations and 
rotations in the plane. In the presence of circulation, new gyroscopic forces enter the equations and 
the Hamiltonian description of the equations involves a non-canonical Poisson bracket on se(2)*. 

In section |3.3[ we define the cyclic centroid of a rigid body with circulation and we show how this 
distinguished point can be described in terms of conformal geometry of the body. Apart from the 
material in the last section and the discussion of the Poisson structures, which is new (to the best of 
our knowledge), most of the material covered in the first two sections can be found, for instance, in 
|26^ [Hj as well as in the classical works of Lamb [T9] and Milne-Thomson j24j . 



3.1 Rigid body motion in potential flow 

The total kinetic energy, T, of the solid-fiuid system (excluding the circulatory motion) is the sum 
of the kinetic energy Tg of the rigid body and the energy Tjr of the fiuid. As both and Tjr are 
quadratic functions on se(2), so is the total energy T. In the absence of external forces or circulation, the 
Lagrangian L of the solid-fluid system is just the kinetic energy: L = T, and the associated equations 
of motion are the so-called Euler-Poincare equations [23] on se(2). These equations determine the 
evolution of the body linear and angular velocity ^ = (w, fi, ^2), and from the reconstruction equations 



(2.5), the position and orientation of the body can then obtained as a function of time. The latter 
determine a curve in SE(2), and it can be shown that this curve is a geodesic with respect to the 
(left-invariant) kinetic energy metric. 

The function L on se(2) is called the reduced Lagrangian. Explicitly, it is given by 

m = \eK, 

where ^ = {oj,vi,V2Y' ^ se(2) = is thought as a column vector and the matrix I is the sum of the 
inertia matrix I_b of the rigid body and the added masses and inertia Ijr of the fluid: I = + Ijr. 
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For future reference, we note that L induces a function C on TSE(2) by left translation, defined by 



^{9:9) = L{9 ^9); see also 



We now focus on the Hamiltonian formulation of the equations of motion for a rigid body in 
potential flow. Recall that a typical element of the dual sc(2)* is written as a row vector fj, = {k,pi,p2) 
and represents the angular and linear momenta of the rigid body (in the body frame). With this 
identification, the Legendre transform associated to L is defined as the mapping 

FL :se(2) ^ se(2)* 

given by FL(^) = /i, where jj. = (I^)^- With the notations of Figure [2| the components of ;U = {k,pi,p2) 
are explicitly given by 

k = {Zb + mia^ + b^) +ljr)uj + {-mb + }Cjr)vi + (ma + ICjr)v2, 
= {-mb + /C^)a; + (m + Mj^)vi + M^^V2, (3.12) 
P2 = (ma + ICjrp + M^'^vi + (m + Mf)v2. 

In classical hydrodynamics k and ipi,P2) are known as "impulsive pair" and "impulsive force" respec- 
tively. The reduced Hamiltonian H : se(2)* — t- M is given by 

H{^,) = ^/xrV^, (3.13) 

and the corresponding (minus) Lie-Poisson equations are /i = ad*„i^/i, where the infinitesimal coadjoint 
action is given by ( |2.7| ) . Written out in component form, these equations are nothing but the Kirchhoff 
equations: 

k = V2Pl-ViP2, ^^^^^ 
Pi=ujp2, p2 = -ujpi, 

where the velocities (w, vi, ^2)"^ and the impulses {k,pi,p2) are related by the Legendre transformation 



(3.12). Finally, we remark that the motion of the body in space can be found from a solution of (3.14) 



by solving the reconstruction equations (2.5) 



3.2 Rigid body motion with circulation 

In the presence of circulation, the Kirchhoff equations on 5e(2)* have to be modified to include the 
Kutta-Zhukowski force. This is a gyroscopic force, which is proportional to the circulation n. In this 
case, the equations of motion are referred to as the Chaply gin- Lamb equations [6l [19] , and they are 
given by 

k = V2P1 - V1P2 - p{avi + l3v2), 
Pi = ^P2 - i^pv2 + pau), (3.15) 
P2 = -ojpi + Kpvi + p(3oj, 
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The constants a and /3 are proportional to the circulation k and depend on the position and orientation 
of the body axes. They are explicitly given by: 



a 



XV(j)Q ■ dl, 



/3 



dB 



YV(j)o ■ dl, 



dB 



(3.16) 



where, as before, {X, Y) are body coordinates for material points in dB. The Chaplygin-Lamb equations 
are discussed in detail in [2j and a derivation from first principles, using techniques from symplectic 
geometry and reduction theory, is given in |35j. 

One easily verifies that if the center of the body axes is displaced to the point with body coordinates 
(r, s), so that the new body coordinates are X = X — r, Y = Y — s, then the circulation constants 
relative to the new coordinate axes take the form a = a — rn, j3 = j3 — sk. Thus, there is a unique 
choice of the body axes that makes these constants vanish. On the other hand, it is also often desirable 
to choose the body axes so that the total inertia tensor I is diagonal. For an asymmetric body, these 
two choices are in general incompatible, see e.g. [19] . 

For our purposes, the choice of body axes will be made to simplify the expression of the nonholo- 



nomic constraint. We therefore consider equations (3.15) in their full generality where a, /3 7^ 0, and 

/? = and 



I is not diagonal. This contrasts with the treatment in [35] where it is assumed that a 
with [2j where the complementary assumption, namely that I is diagonal, is made. 



A calculation shows that equations (3.15) are Hamiltonian with respect to the usual Hamiltonian 
H : se(2)* — )■ R given in (3.13) and with respect to the following bracket of functions on se(2)*: 



{F,i^},,„,/3(/u) = (V^F) 






iiV.F) X (V^K)) . 



(3.17) 



where, "x" denotes the standard vector product in M'^, and, as before, V^F is the gradient of F with 
respect to the variables {k,pi,p2) = /i. In principle one would need to show that the above bracket 
satisfies the Jacobi identity. This will follow from the results of section [4] where it will be shown that 
the above bracket is Lie-Poisson on the dual Lie algebra of a central extension of SE(2) by M.^. 



The first term in (3.17) is the usual (minus) Lie-Poisson bracket on se(2)* and defines the dynamics 
in the absence of circulation. The second term, which vanishes in the absence of circulation, is the 
remnant of the magnetic term in the reduced form in T* SE(2). For future reference, we record here 
the form of the bracket for a = f3 = 0: 



{F,K},ifj.) = iV,F) 





dF dK 
dpi dp2 



dF dK 
dp2 dpi 



(3.18) 



3.3 Conformal geometry and the constants a,/3 

In this section, we describe some of the complex geometry underlying the constants a and f3 defined in 
(3.16). As discussed in the previous section, there exists a unique choice of origin for the body frame 
relative to which a and (3 vanish. We denote this point by Tc = {xc, Vc) and refer to it as the cyclic 
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centroid. Using (2.9) and (3.16) one deduces 



1 



K Jqq dl 

In general, the location of this point will be different from the geometric centroid of the body, but 
it shares one important characteristic with the centroid: it can be found among the fixed points of the 
isometrics of the body, if any. As an easy consequence, we have that for a homogeneous rigid body with 
two symmetry axes, the cyclic centroid, the geometric centroid, and the center of mass all coincide. 
For asymmetric bodies, the three points will in general be different. 

Proposition 3.1. Let T : B ^ B he an isometry of the body, i.e. a linear isomorphism o/ such 
that T{B) = B. Then Yc is a fixed point ofT. 



Proof. Using the linearity of T we have 



T{rc) = -f T{r)—dl = -(b r —dl = - f r —dl , 

K Jgjg dl K JT-UdB) 91' K /^B 81' 



since T leaves dB invariant, and therefore T{rc) = Tc. □ 

The location of the cyclic centroid can be made more explicit using some methods from complex 
geometry. Let W{z) = (po + iipo be the complex potential associated to the cyclic motion and put 
Zc = Xc + iyc- We then have that along the boundary 

-^dl = 3? I — — dz] and -^v-dl = 9 ( — — dz 
dl \ dz J ol \ dz 

but since the boundary of the body is a streamline, the imaginary part of ^^(z) is constant on the 
boundary, and hence 

1 f ^ (dW . \ 1 



^ z^[—dz]=-(p zdW, 
1^ Job \ dz J K J 

where the integral on the right-hand side is taken along an arbitrary contour around the body. 

Now choose an arbitrary conformal map ( z = z{Q from the interior of the unit disc to the 
exterior of the rigid body, and assume that z{0) = oo. The complex potential is related to the 

complex potential tf(C) in the ("-domain by VF(z) = w{C,), and moreover w(C) = 2^1o§C- As a result, 
the cyclic centroid can be expressed as a contour integral along the unit circle as 



1^ /|C|=i 27ri C 

Evaluating this contour integral requires us to know the constant term oq in the Laurent series for 
around 0. Since only has a simple pole at infinity, we have that Zc = —oq. 

These manipulations hint at a neat geometric interpretation for Zc in terms of the first Faber 
polynomial and the capacity of the rigid body, viewed as a domain in the complex plane |8l|9]. Recall 
that the Faber polynomials of a complex domain D arise when considering a uniformizing map $ : z i— )■ 
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Figure 4: Location of the cyclic center (red dot) for a triangular and an irregular body. This point 
is the center of the level set |/i(z)| = 1 (circle indicated in green) of the first Faber polynomial 
fi{z). The radius of this circle is the conformal capacity of the body. 



w := ^{z) from the exterior of D to the exterior of the unit disc. After suitable normalizations, this 
map has a Laurent expansion 

^{z) = Y^{z + ao + aiz~^ + a2-z~^ H ) 

1^1 

where C is the capacity of D. The mth Faber polynomial fm{z) is then the polynomial part of [<l>(z)]'". 
In particular, we have that fi{z) = \C\~^{z + oq). 

Now, note that the level set |^*(^;)| = 1 is precisely the boundary of the rigid body. In the theory 
of complex approximation, this level set is approximated by the level sets |/m(-2)| = 1; which provide 
successively better approximations as m increases. The first nontrivial level set, |/i(-z)| = 1, is a circle 
whose center is precisely the cyclic centroid Zc and its radius is the capacity \C\ of the domain. From 
this point of view, the cyclic centroid is the center of the best-fitting circle through the boundary of 
the rigid body. In Figure [4| a few examples are shown of polygonal rigid bodies, for which z{C,) is a 
Schwartz-Christoffel mapping (see |9j). Also indicated is the level set of |/i(z)| = 1. 



4 Geodesic flow on central extensions of Lie groups 

In this section, we give a geometric interpretation of the Chaplygin-Lamb equations. In particular. 



we show that the non-canonical bracket (3.17) on se(2)* can be made into a canonical bracket on a 
"larger" space, obtained by considering a central extension of SE(2). After reviewing the geometry of 
central extensions, we relate the particular extension, and the cocycles used in its description, to the 
circulation constants k, a, and (3. 



After constructing the relevant central extension G of SE(2), we show that the bracket (3.17) is 
precisely the Lie-Poisson bracket on the dual g* of the Lie algebra, and we argue that the Chaplygin- 
Lamb equations hence give rise to geodesic flow on G. This result is interesting in its own right, and 
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provides the first step towards showing that the hydrodynamic Chaplygin sleigh with circulation can 
be formulated as an LL-systeni, i.e. a nonholonomic system on G whose Lagrangian and constraints 
are invariant under the left action of G on itself. 

To construct the central extension G, we proceed in two steps. We first assume that a = (3 = 
(recall that this can be achieved by a suitable choice of origin for the body- fixed frame), and we recall 
from [35] that the circulation k is described by means of a real- valued cocycle Ci, giving rise to a 
central extension of SE(2) by M, known as the oscillator group [33\. The effect of nonzero a and (3 can 
then be described by a second cocycle C2 with values in M?, giving rise to a central extension of the 
oscillator group by M^. By combining both of these cocycles, we then describe G as a central extension 
of SE(2) by M3. 

In terms of Lie algebra cohomology, the cocycle Ci is closed but not exact, while C2 will be shown 
to be exact. It is known that central extensions by exact cocycles are trivial, in the sense that the effect 
of C2 can be "gauged away" by a suitable redefinition of the momenta of the system (this is equivalent 
to the fact that a and /3 can be made zero by a suitable choice of origin), while the circulation k 
(described in terms of the non-exact cocycle Ci) cannot be made to vanish. 

4.1 Central extensions of Lie groups 

We briefly recall the definitions and basic properties of central extensions of Lie groups to introduce the 
relevant notation. A detailed account of the geometry of central extensions in the context of mechanics 
can be found in [22l EH [l6] . 

Definition. Let G be a Lie group and A an abelian Lie group. We will use additive notation for the 
group operation in A. Let G be an extension of G by A, that is, there is an exact sequence 

0^ A^ G ^ G ^ {e}. 

Throughout the rest of this paper, we will assume that G = G x ^ as a manifold and that the group 
multiplication on G is 

{g, a){h, /3) = {gh, a + P + B{g, h)), (4.19) 

where B : G x G ^ A is a. normalized group two-cocycle. Associativity of the multiplication on G is 
equivalent to the two-cocycle identity, 

B{f,g) + B{fg,h) = B{f,gh)+B{g,h) for aU f,g,heG. 

The assumption that the two-cocycle B is normalized can be made without loss of generality and 
amounts to 

B{g, e) = B{e, g) = for aU g e G, 

implying that B{g,g^^) = B{g^^,g). Under these assumptions, G is referred to as a central extension 
ofG by A. 
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Central extension of Lie algebras. The Lie algebra g of G is isomorpliic as a vector space to g x a, 
and is equipped with the following bracket: 



[(C,a),(r/,%=(K,^]g,C(6r?)) 
for a), (ry, b) G Q. Here, the a- valued Lie algebra two-cocycle C : g x g 



dsdt 



{B{g{t),h{s))-B{h{s),gm, 



a is defined by 
i = l,2, 



s=t=0 



with g{t), h{s) smooth curves on G satisfying g{0) = ^, h{0) = rj. 

It may happen that C can be written in terms of a one-cocycle A : g 
for all S,,Tj £ Q. In this case, C is a coboundary in the sense of Lie algebra cohomology, and we write 
C = 6A. When this happens, the central extension is said to be trivial, as the mapping 



a as C(^, rj) 



(4.20) 



^^:gxo^g, {C,a) ^ {C,a - A{C)) 



(4.21) 



then determines a Lie algebra isomorphism between the Lie algebra g x a with the product bracket and 
the central extension g. 



Lie-Poisson structures. As a vector space, the dual Lie algebra g* equals g* x a*. For (^,0") E g*, 
the lb Lie-Poisson bracket of functions F,K£ (g* ) is readily computed to be 



{F,i^}±(^,a) 



6F 6K 
5^ ' Sfj, 



(4.22) 



Notice that this bracket only involves functional derivatives with respect to /x. Therefore, if we fix 



the value of o" G a* in (4.22), we obtain a non-canonical Poisson bracket on g*, given by formally the 



same expression as (4.22): for /, A; G C°°(g*) we have 



{f,k}tip) = ±(fi 



5f 6k 



6fi^ 6fj, J 



(4.23) 



where o" G a* is now regarded as fixed. Roughly speaking, we therefore have a one-to-one correspondence 
between Poisson brackets on g* which are the sum of a Lie-Poisson term and a cocycle, and Lie- 
Poisson brackets on central extensions g*. This can be made rigorous by observing that the injection 
io- : g* ^ g*, given by io-(//) := (/^, c) for a fixed, is a Poisson map taking the non-canonical Poisson 



structure (4.23) into the Lie-Poisson structure (4.22) 



In the case of a trivial central extension, the cocycle term in (4.23) can be "gauged away": let 
= 6A. For a G a*, we denote by Aa : g — )• M the linear map defined by Aa{(,) ■= {cr,A{£,)) for all 



C 

^ G g. Note that Aa G g*. The Poisson bracket (4.22) can then be rewritten as 



{F,K}f,i^,a) 



±(fl-Aa, 



t((^,A 



6F 6K 
6fi ' 6fj, 



'5F 6K' 

'6F 6K 
6^ ' 5fi 

so that the cocycle vanishes apart from a shift ^ — Aa in the momenta of the system. 
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Proposition 4.1. The shift map : g* — )• g* x a*, given by $_4(/x, a) = {fi — Aa,cr) is the dual of 



the trivialization mapping (4.21). Moreover, the shift map is a Poisson map, taking the Lie-Poisson 



structure (4.22) on g* into the product Lie-Poisson structure on g* x a* 



Hamiltonian vector fields. The Hamiltonian vector field of a function H on g* is defined by the 
equations 

/i = =F ^adlff/i + (7 o C -^^ , o- = 0. 

Notice that the components of a are Casimirs of the bracket and as such they are constants of motion, 
in our case these constants of motion will be (a multiple of) the circulation k and the constants a, (3. 



4.2 The oscillator group 

We start by defining the real valued SE(2)-two-cocycle Bi : SE(2) x SE(2) 



by 



(4.24) 



where J is the 2x2 symplectic matrix 



1 
-1 



This cocycle differs from the one used in |35j by a multiplicative factor of ^. On the Lie algebra level, 
we have by (4.20) that the infinitesimal cocycle Ci : se(2) x se(2) — t- M is given by 

Cl ((W, Vl,V2), (w', v[,V2)) = Viv'2 - V2V[. 

The central extension of SE(2) by M using the cocycle Bi is referred to as the oscillator group [33]. 
and will be denoted as Osc. The Lie algebra 05C of the oscillator group is isomorphic to se(2) x M, with 
Lie bracket 

[{(jJ,Vl,V2;z) , {J ,v'i,v'2;z')]^^^ = {\^{uj,Vi,V2) , (w',fi,f2)]^^(2) ; '^1 {{^,Vl,V2),{uj' ,v[,v'2))^ 

= ( , V2ijj' — iov'2 , ^^v'y — VlUj' ; viv'2 — V2v'i ). 



On the dual osc* = se(2)* x M*, the (minus) Lie-Poisson bracket is given by (4.22), or explicitly by 



{F,i^}„-,.(/i,ao) = - ^oCi ^) 



for (//, (To) G 5e(2)* x M*. It is easy to see that this bracket is precisely the bracket (3.18), where ctq 
plays the role of the circulation pK. We conclude that the effect of non-zero circulation can be described 



in terms of the geometry of the oscillator group, and more precisely the cocycle (4.24). 



In terms of Lie algebra cohomology, it is easy to show that Ci is a closed two-cocycle which is not 
exact. Furthermore, as //^(se(2),M) is isomorphic to M, we have that Ci determines a generator of 
the second cohomology. Since isomorphy classes of central extensions are classified by the second Lie 
algebra cohomology, we will refer to the oscillator group as the central extension of SE(2) by M. 
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4.3 A central extension of SE(2) by 



We now describe the effect of nonzero a and /3 in the Poisson bracket (3.17). To this end, we observe 
that the equations of motion (3.15) for nonzero q and /3 can be obtained from the equations where 

(4.25) 



a = /3 = by making the substitution, or momentum shift, 



where cJi = — p/3 and o"2 = pa, while the angular momentum k remains invariant. This momentum shift 
can be described in geometric terms as follows. Let A : se(2) — t- be the linear map A{uj, vi,V2) = 
(0, vi, V2)- This map is a one-cocycle on se(2) with values in M'^ and its derivative is given by 

SA{{u;,vi,V2),{uj',v[,V2)) = -A{[{uj,vi,V2), {uj' ,v[,V2)]s,{2)) = {0 , UJV2 - V2Uj' , viUj' -u;v[). 
Now, consider fixed values cri,o"2 G M and define the linear map ^(0-1.0-2) • ^^(2) ~^ ^ given by 



A 



(0-1,0-2 



){u},vi,V2) := ((0, cJi, 0-2), ^(a;,f 1,^2)) = criVi + cr2V2. 



As ^(0-1,0-2) is an element of se(2)*, we may write the momentum shift (4.25) more formally as the map 
: se(2)* se(2)* given by 

^.4(/^) = /^- Aoi,02), (4.26) 
for all fj, G se(2)*. The minus sign is due to the fact that this is the active version of the transformation 
(4.25): if we denote the new momenta of the system by pi,p2, then the effect of performing the 
substitution (4.25) is that the old and the new momenta are related by pi = pi + (7i,p2 = P2 + cr2, 
which is just (4.26). 

As the dual of the Lie algebra of the oscillator group is just the Cartesian product se(2)* x M*, the 
map <I>yi gives rise to a map (/i, ctq) i— )• (^'yi(/^), o"o) on osc* which we denote by as well. We now 
investigate the behavior of the Poisson structure (3.18) under the map <I>^. Since is a constant 
shift map, we have for arbitrary functions F, K on osc* that 

5{F o _ 5F 

and similarly for K, so that {F o ^^^,K o {p, cr) = {F, K}'^^^, (/x, a). On the other hand. 



{F,K}-^,{^^\p),a) = {F,K}-,^.,i<^>^\f,)) - aoCi 



«(2)* 
{^'^}«(2).(/^) 



A 



(01,02); 



^\ 
\6fi' 6fiJ 
'6F 

6fi ' 6n 



o"oCi 



6F 6K 



( 



6F 6K 



(4.27) 



so that $_4 takes the Lie- Poisson bracket on osc* into the bracket (4.27) modified by a cocycle. The 
cocycle aoCi — 0-2) can be made more explicit by noting that 



(4.28) 









V1V2 


- V2v'i 








'J' 




(cJoCi - 5A(^„^^„^)){{u;,vi,V2), {uj',v[,V2)) = 






V2<jj' 


- ^V'2 


-A 


Vl 


X 




) 




^2. 




ujv'i 


— VlOj' 




.^2. 




A. 
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where cr := (cto, (Xi, (T2)- Note that this is precisely the cocycle in the Poisson bracket (3.17) for the 
specific value of <t := {pK, —pj3,pa). We conclude that the eff'ect of introducing nonzero ai,(T2 is to 
modify the Lie- Poisson bracket on osc* by a trivial M^-valued cocycle —5A(^„^^„^y We now let C2 be 
equal to —5 A: this is an M'^-valued cocycle on se(2), which can be integrated to a group two-cocycle 
B2 on SE(2), given by 

S2((i?e,x),(i?esx')) = (0,i?ex'-x'). (4.29) 



As we showed previously, we can describe the cocycle bracket ( |4.27 ) by considering the central 
extension of the oscillator group by the cocycle B2 given in (4.29). This is equivalent to extending 
SE(2) by using the combined cocycle B := i?2) : SE(2) x SE(2) — ^ M^; the result is an extension 
G which is isomorphic to SE(2) x with multiplication 

(i?6i,x; w) • {R0i,yi;W) = {Rg+0/,x + figx ; + w' + B((i?e,x), {Rg>,x'))). 

The infinitesimal cocycle C : 5e(2) x 5e(2) — t- associated to B is given by C = (Ci, C2), or 

C{{U}, Vl,V2), {uj', v'i,V2)) = {VIV2 — V2v'i,V2i^' — a;w2, ww'i — viu') 



{UJ, Vl,V2) X {uj',v[,v'2), 



(4.30) 



which is precisely the cocycle appearing on the right-hand side of ( |4.28 ). The bracket on the algebra 
Q := Lie(G) is then given by 

[{uj,vi,V2;z),{uj',v[,V2;z')]'^= (^[{u},vi,V2) , (a;',ui,U2)]s,(2) 5 C {i^,vi,V2), {uj' ,v[,v!2))^ 

= ( , V2U^' — LOv'2 , U^v'i — Vlbj' ; Viv'2 — V2v'i,V2Uj' — UJV'2,UJV[ — VlLo' ) . 



The Lie-Poisson bracket on q*. We write explicitly the (minus) Lie-Poisson bracket on the dual 
Lie algebra g*. Note first that as a vector space g* is just se(2)* x M^, so that an element u of g* can 



be written as u := (^, cr) where p, S -Se(2)* and cr E M^. In view of (4.22) and (4.30) we obtain: 
Proposition 4.2. The (minus) Lie-Poisson bracket of functions F,K£ C°°(q*) is given by 



(5F 5K\ 



{F,K}-(p,a) = {F,K}-,^Ap)-a 



where we identify se(2) with in the usual way to make sense of the vector product of the functional 
derivatives. 

In coordinates, the Poisson structure of the above proposition is 



{F,i^K(/x,cT) = (V(^,^)F) 



T 








-P2 - 0-2 


Pl + 0"! 








\ 




P2 + 0-2 





-do 













-Pl - 0"! 
























































v 

















/ 



(V(^,,.)i^), 



where as usual p = {k,pi,p2) G se(2)* and cr = (ctq, cji, 172) G 
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The equations of motion. The reduced kinetic energy Hamiltonian H{n) = ^fi ■ I ^/x on se(2)* 
defined previously in (3.13) gives rise to a Hamiltonian H^* on g* as follows: 



H^,ils,a)=Hif,) + -\\(7\ 



(4.31) 



The equations of motion for this Hamiltonian and the bracket {•, ■}^, are then given in coordinates by 



k = —p2Vi +P1V2 — (J2V1 + aiV2, 
Pi = + 0-20; - aoV2, 
P2 = -Pi^ - ci"^ + o'ot'i, 

d- = 0, 



that coincide with (3.15) if we put o"o = /OK, ui = —pP, (T2 = pct- We mention that, in addition to the 



energy H, these equations possess the conserved quantity 

F(/i, cr) =pI+pI + 2aok + 2aipi + 2a2P2, 

which is a Casimir function of the bracket {•, ■}^, . In fact, for a constant value of cr, the regular level sets 
of F are symplectic leaves of the symplectic foliation of g* . It is easily seen that these are paraboloids 
if ctq 7^ 0, and cylinders otherwise. The trajectories of the system are defined by the intersection of the 
level sets of F and H. 



Notice that the Hamiltonian (4.31) is the (reduced) Legendre transform of the left G-invariant 
kinetic energy Lagrangian Cq : TG — M whose value at the identity is Lg : g — >• M given by 



^(e,a) 



1 1,, , 

+ - a 



for (,^,a) e se(2) x 



Our discussion can be summarized in the following theorem. 



Theorem 4.3. The Chaplygin-Lamb equations (3.15) are of Lie-Poisson type on the dual Lie algebra 
g* of the central extension G o/SE(2) byM? defined above, and with respect to the left invariant, kinetic 



energy Hamiltonian determined by (4.31). In particular, the bracket (3.17) coincides with the restriction 
of the Lie-Poisson bracket on g* to the Casimir level set defined by ctq = pi^-, = ~pl^-, ^2 = pct. 



5 Nonholonomic LL systems on central extensions 

We have seen that the Chaplygin-Lamb equations for the motion of a rigid body in a potential fluid in 
the presence of circulation are Lie-Poisson equations on g* corresponding to a pure kinetic energy left 
invariant Lagrangian. Our goal in this section is to add a left invariant nonholonomic constraint to the 
system. The resulting reduced equations, which are consistent with Lagrange-D'Alembert's principle 
stating that the constraint force performs no work during the motion, are the so-called Euler-Poincare- 
Suslov (EPS) equations. 
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5.1 The case of a general central extension 



To the best of our knowledge, the derivation of the EPS equations has never been made exphcit in the 
case where the underlying Lie group is a central extension so we take the opportunity to develop the 



general theory in detail. In section 5.2, we apply this theory to the hydrodynamic Chaplygin sleigh. 



Euler-Poincare-Suslov equations on a Lie group G. In general, a nonholonomic system on a 
Lie group G with a left invariant kinetic energy Lagrangian and left invariant constraints is termed an 
LL system. Due to invariance, the dynamics reduce to the Lie algebra g, or to its dual q* if working 
with the momentum formulation. We start from a reduced Lagrangian L : g — J- R, which defines an 
inertia operator I : g — )• g* by the relation 

L(e) = ^(ie,o, forees, 

where (•, •) denotes the duality pairing. The reduced Hamiltonian, : g* — )■ M, is then given by 

^(/^) = ^(/^>I~V>, for ^ eg*. 

The nonholonomic constraints are expressed in terms of n linearly independent fixed covectors 
G g*, i = 1, . . . , n: we say that an instantaneous velocity ^ G g satisfies the constraints if 



= 0, 



1, . . . ,n. 



(5.32) 



We let D C g be the vector subspace of all velocities satisfying the constraints, and we say that the 
constraints are nonholonomic if T) is not a Lie subalgebra of g. 

The reduced EPS equations on g* are given by, see e.g. |3], 



fjL = adj_i^/i + ^Aji/i, 



(5.33) 



i=l 



where the multipliers Aj, i = 1, . . . , n, are certain scalars that are uniquely determined by the condition 



that the constraints (5.32) are satisfied. Explicitly, the Lagrange multipliers are given by 

n 



A. 



(5.34) 



i=l 



where D is the matrix with components Dij := ^i^j), hi = 1, . . .n, and D ^ is its inverse. 



Euler-Poincare-Suslov equations on central extensions. Now suppose that G is a central ex- 
tension of G by the abelian Lie group A as explained in section kl As before, we assume that G = G x A 



with the multiplication (4.19). We let L be a left-invariant Lagrangian on G, with associated Hamil- 
tonian H, and we let Vi £ g*, i = 1, ... ,n be a set of linearly independent constraint covectors. We 
now wish to "lift" these data to the central extension G, so that we can derive the corresponding EPS 
equations on the co-algebra g* . 
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By left translating the co- vectors Vi G g* one can define left-invariant constraint one- forms e^, on G, 
given by ei{g) = T*Lg-iVi £ T*G, i = 1, . . . ,n. Since G = G x A, these constraint one-forms naturally 
induce constraint one-forms ej on G, given by ti{g,a) = {e{g),0). Likewise, the co-vectors Ui G g* can 
be lifted to co- vectors z/j = (fj, 0) in g*, and we have that 'ei{g, a) = T^^ ^^L(g cf)-iz7j. 

Secondly, we define the left invariant, kinetic energy Hamiltonian Hrp,Q : r*G — )• M, whose value 
at the identity is given by 



where || • |p denotes any positive definite, quadratic form on a*. For convenience we write 

i7g,(/i,cj) = ( (^,cj) , r\^,cr) 



(5.35) 



where the non-degenerate, extended, inertia tensor I : g — )• g* is determined from (5.35). 

We have thus extended the left invariant constraints and kinetic energy on G to define an LL system 
on the central extension G. The corresponding EPS equations (5.33) become, 

n 

ifi, &) = adYi(^^^)(/i, cr) + Y^ 

i=l 



i=l 

where ad* (respectively, ad*) denotes the infinitesimal coadjoint action on g* (respectively, on g*). 
Here, the Lie algebra two-cocycle C : g x g — )• a is defined by (4.20). 

The above equations can be written as 



ad*_i^;U + Cr O G(I V: ■ ) + ^i'^i^ CT = 0, 



(5.36) 



i=l 



and one finds the following expression for the multipliers: 



Y,iD'%j ((ad*_i^M,r'^i> + (c^,c(rV, r^i^O) 



(5.37) 



i=l 



Notice that these equations reduce to (5.33), (5.34) if a = 0. 



Existence of invariant measures. It is natural to ask whether the equations (5.36) and (5.37) 
possess an invariant measure. We settle this question using the criterion of Jovanovic [Ej for the 
existence of an invariant measure for the EPS equations on the Lie algebra of a non-compact group. 
We assume that there is only one constraint, so that n = 1; the case of multiple constraints can be 
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dealt with in a similar way. Following [T^, the necessary and sufficient condition for equations (|5.32|), 

at the constraint co 

ad*j_i^z? + f = cv, 



(5.33) to have an invariant measure is that the constraint covector v = ui G 3* satisfies 

1 



c G 



(5.38) 



where T G g* is defined by the relation (T, = trace(ad^), ^ G g. 
However, for a), (r/, 6) G g x a = g we have 

ad(g,a) = (ad^?/,C'(?,??))- 
It follows that the operator ad(^ a) :g = gxa— ;'g = gxa has matrix representation 



ad^ 



Hence trace(ad(^ a)) = trace(ad^) and we can write T = (T, 0) G g* x a*, where T G g* is defined by 
(T, ^) = trace(ad^) for ^ G g. In addition, since v = (i^, 0), we can write (5.38) in components as 



-^^j^p^ad*(i-ii,^o)(i^,0) + (r,0) = c(z^,0). 



c G 



Therefore, the condition (5.38) is equivalent to 

/ TT-i \ ad*i-i^t^ + T 



c G 



which is precisely the necessary and sufficient condition for the equations (5.33), (5.34) to possess an 
invariant measure. This analysis can be generalized to the case where the number n of constraints is 
arbitrary. This shows: 



Theorem 5.1. The Euler-Poincare-Suslov equations (5.36) and (5.37) on the dual Lie algebra g of 
the central extension G of the Lie group G, possess an invariant measure for an arbitrary value of 
a E a* if and only if they possess an invariant measure for the specific value a = 0. In other words, 
such an invariant measure exists if and only if the Euler-Poincare-Suslov equations (5.33), (5.34) on 
g* possess an invariant measure. 



5.2 The hydro dynamic Chaplygin sleigh with circulation 

We are now ready to consider the mechanical system of our interest which is the generalization of 
the hydrodynamic version of the Chaplygin sleigh treated in [TO] to the case when there is circula- 
tion around the body. Recall that the classical Chaplygin sleigh problem (going back to 1911, [7]) 
describes the motion of a planar rigid body with a knife edge (a blade) sliding on a horizontal plane. 
The nonholonomic constraint forbids the motion in the direction perpendicular to the blade. In its 
hydrodynamic version, the body is surrounded by a potential fiuid and the nonholonomic constraint 
models the effect of a very effective fin or keel, see |10J. 
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With the notation from section [2| we let {Ei,E2} be a body frame located at the contact point of 
the sleigh and the plane, and so that the Ei-axis is aligned with the blade (see Figure [T]). The resulting 
nonholonomic constraint is given by V2 = 0, and is clearly left invariant under the action of SE(2), as 
it is solely written in terms of the velocity of the body as seen in the body frame. 

In the absence of constraints, the motion of the body is described by the Chaplygin-Lamb equations 



(3.15) that, as we have shown (Theorem 4.3), are Lie-Poisson equations on the dual Lie algebra of the 
central extension G of SE(2) and with respect to a left invariant kinetic energy Hamiltonian. The 
reduced nonholonomic equations are thus of EPS type on q* . Note that the co- vector (0,0,1) G q* 
annihilates all elements ^ = {uj, vi,V2) G se(2) that satisfy the constraints. Thus, in agreement with 
the results of Theorem 4.3, by putting /i = {k,pi,p2) £ se(2)* and a = (pn, —pf3, pa) E M^, we get the 



following explicit expression for the reduced nonholonomic equations (5.36), 



k = V2P1 - V1P2 - p{avi + l3v2), 
pi = ujp2 - Kpv2 + pau}, 
P2 = -ojpi + Kpvi + p(3uj + A, 

where the multiplier A is determined from the condition ^2 = 0. 



(5.39) 



The total inertia tensor and the circulation constants a, /?. The expression for Ig with respect 



to the body frame {Ei, E2} is given by (2.6) where m is the mass of the body, (a, b) are body coordinates 
of the center of mass, and I is the moment of inertia of the body about the center of mass. While 
this simple expression is independent of the body shape, an explicit expression for the tensor of adjoint 
masses Ijr can be given explicitly only for rather simple geometries. A simple yet interesting one is 
that for an elliptical uniform planar body with the semi-axes of length A > B > 0. Assume that the 
origin has coordinates (r, s) with respect to the frame that is aligned with the principal axes of the 
ellipse and that the coordinate axes Ei E2 are not aligned with the axes of the ellipse, forming an angle 
V (measured counter-clockwise) , as illustrated in Figure [Tj 

For this geometry, starting from the formula ( 2.11[ ) for the added inertia tensor in a body frame 
whose axes are aligned with the symmetry axes of the ellipsoid, one can show that 



Ijr = p-K 



I 



\ 



■^^ 2 — + s { B COS \i ^ A sm ti) 

+ r^(A^ cos^ u + sin^ u) - rs(A^ - S^) 3111(2^) 

-lr{A^ - S2)sin(2u) 

— r{A? CO!? V + sii? v) 
+ 1s(A2 - B^) sin{2ij) 



s{B^ cos^ V + A'^ sii? -u) 
-^r(A^ - B^) sin{2ij) 

cos^ V + A^ sii? D 

A^-b'^ ■ ,„ , 
2 sin(2^j) 



-r(A^ cos^ V + B^ sin^ v) \ 
+ 5s{A^ - S^) sin{2u) 



n{2ij) 



A^ cos^ i> + B^ sin^ v 



. (5.40) 



The total inertia tensor, I = + Ij-, of the fluid-body system is then given by 



/ I + m(a^ + b^) 

+P-!tI '-^ ' + s'^(B^ cos^ v + A^ sin^ u) 

+ (A^ CO? D + sin^ v) - rs(A'^ - S^) sin(2u)] 



— mb + p'!t[s(B^ cos*^ v + A'^ sin"^ ti) 
-^r{A^ - B^) sin(2-u)l 

ma + p7r[ — r(A^ cos^ v + B^ sin"^ ti) 
+ ls(A^ - B^)Bin(2v)] 



-mb + p7t[s{B cos + A sin u) 



m + p7r[B'^ cos^ v + sin^ u] 
pvr I J sin(2i>)l 



ma + p7r[ — r{A cos v B sin v) 
+ iB(A^ - B^) sin(2-u)l 



(A^-B^ . ,„ ,1 
pvr 5 sin(2ij) 



m + p7v(A^ CQs'^ V + B^ su? v) 
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Notice that in the presence of the fluid, if f 7^ n|, n G Z, the coefficient I23 = I32 is non-zero. This 
can never be the case if the sleigh is moving in vacuum as one can see from the expression given for 



in (2.6). The appearance of this non-zero term leads to interesting dynamics that, in the absence of 



circulation were studied in [T^. For the above geometry, the circulation constants a and /3 defined by 



(3.16) are computed to be: 

a = K,{rcosv + ssinu), /3 = k(— r sinf -|- scosv). 

Notice that, in the presence of circulation, the two constants, a and /3, can vanish simultaneously only 
if r = s = 0, that is, only if the contact point is at the center of the ellipse. 

In the sequel we assume that the shape of the sleigh is arbitrary convex and that its center of mass 
does not necessarily coincide with the origin, which leads to the general total inertia tensor 



J 


-L2 


Li 


-L2 


M 


Z 




Z 


N 



and with arbitrary circulation constants q, /3. 



Detailed equations of motion. In accordance with ( |5.34 ), the multiplier A is given by: 

1 



A 




V2P1 — V1P2 — p{avi + fiv2) 
ujp2 — Kpv2 + pauj 
-ujpi + npvi + pj3ijj) 



where 



I" 



1 



det(I) 



MN - ^2 ZLi + NL2 
ZLi + NL2 JN - Lf 
-ZL2 - MLi -L1L2 - JZ 



-ZL2 - MLi 
-L1L2 - JZ 
JM -Ll 



A long but straightforward calculation shows that, by expressing uj,vi and V2 in terms of k,pi,p2, 

the constraint V2 = 0, one obtain 

{Lioj + Zvi + pa) {L2U} — Mvi) , 



substituting into (5.39), and enforcing the constraint V2 = 0, one obtains: 

1 



UJ 



Vl 



D 
1 

D 



{Lioj + Zvi + pa) {Joj - L2V1) 



(5.41) 



where we set D = det(I)(I-^)33 = MJ - Lj. Note that D > since I and I ^ are positive definite. 
Note as well that if a = we recover the system with zero circulation treated in so from now on 
we assume a 7^ 0. 

The full motion of the sleigh on the plane is determined by the reconstruction equations which, in 
our case with V2 = 0, reduce to 



X = Vl COS a, 



y = Vl sm ( 
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Figure 5: Reduced phase portrait. The line £ consists of equilibria that are either stable (filled dots) 
or unstable (empty dots). The trajectories are contained in the level surfaces of the Hamiltonian 
H that are ellipses on the Vi u> plane. 



The reduced energy integral has 



1 



H =^ ( Jw^ + Mvl - , 

and its level sets are ellipses on the (w t'i)-plane. 

As seen from the equations, the straight line I = {Liuj + Zvi + pa = 0} consists of equilibrium 
points for the system. Each of these equilibria corresponds to a uniform circular motion on the plane 
along a circumference of radius I — I . 

Notice that if Z = Li = the line ^ of equilibra disappears. In fact, it is shown in [TU] that in 
the absence of circulation the system possesses an invariant measure only for this specific value of the 



parameters. In view of Theorem 5.1 we conclude 



Proposition 5.2. The equations of motion (5.41) possess an invariant measure if and only if Z = 
Li = 0. 

In this case we obtain simple harmonic motion on the reduced plane uj,vi. The reduced phase space 
in the case where Z ^ and Li / is illustrated in Figure [5| 

For any value of the parameters the reduced system ( |5.41[ ) can be checked to be Hamiltonian with 
respect to the following Poisson bracket of functions of w, vi that can be obtained using the results of 

m 

The invariant symplectic leaves consist of the semi-planes separated by the equilibria line i and the 
zero-dimensional leaves formed by the points on this line. 

The integration of these equations as well as a detailed study of the motion of the sleigh on the 
plane will be given in |llj . 
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Conclusions and further work 



In this paper, we presented an interesting example of a nonholonomic hydrodynamical system. We have 
shown how the Hamiltonian structure of the Chaplygin-Lamb equations can be understood in terms 
of central extensions of Lie groups, and used this geometric construction to formulate the equations 
of motion for our problem in EPS form. This structure of the equations was exploited to study the 
existence of invariant measures for the hydrodynamic Chaplygin sleigh with circulation. 

The physical relevance of our model lies in the fact that in the design of underwater vehicles the 
nonholonomic constraint can be interpreted as a simple model for a fin. From the mathematical point 
of view, furthermore, our example is remarkable since its equations of motion are explicitly integrable, 
and describe interesting dynamics as we show in 

For the future, we intend to study the motion of the hydrodynamic Chaplygin sleigh coupled to 
point vortices in the fluid [25] . The equations of motion for interacting point vortices and rigid bodies 
(without nonholonomic constraints) were recently derived in |3H [T] and since then there have been 
significant efforts towards discerning integrability and chaoticity |27| [29] and towards uncovering the 
underlying geometry of these models |34] . We plan on coupling the nonholonomic Chaplygin sleigh 
with one or several point vortices in the flow, taking these models as our starting point. We hope to 
report with progress on these problems in the near future. 
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